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ABSTRACT 

Using high-quality spectra of the twin stars in the XO-2 binary system, we have detected significant 
differences in the chemical composition of their photospheres. The differences correlate strongly with 
the elements’ dust condensation temperature. In XO-2N, volatiles are enhanced by about 0.015 dex 
and refractories are overabundant by up to 0.090 dex. On average, our error bar in relative abundance 
is 0.012 dex. We present an early metal-depletion scenario in which the formation of the gas giant 
planets known to exist around these stars is responsible for a 0.015 dex offset in the abundances of all 
elements while 20 Mq of non-detected rocky objects that formed around XO-2S explain the additional 
refractory-element difference. An alternative explanation involves the late accretion of at least 20 
of planet-like material by XO-2N, allegedly as a result of the migration of the hot Jupiter detected 
around that star. Dust cleansing by a nearby hot star as well as age or Galactic birthplace effects can 
be ruled out as valid explanations for this phenomenon. 

Subject headings: stars: abundances — stars: fundamental parameters — stars: planetary systems — 
stars: individual (XO-2) 


1. INTRODUCTION 

Over the past few years our team has pioneered and 
further developed a technique for measuring elemental 


(e-g-:, 

Melendez et al. 2009 

2012 Ramfrez et al. |2009 

201lj 

2014a|bl 

Bedell et al. 

2014 Liu et al. 2014| I'iticci 

Maia 

et al. 

|20] 

While most other chemical analyses 


u.uu 

-IT 

I o 


we have been 
1 dex (~2 %) 


are limited 

in quantities such as [Fe/H], [0/Fe], etc. 
able to derive relative abundances with 
precision, and even better, down to the 1 % level, in a 
few critical cases. 

The success of our technique can be traced back to two 
key requirements: 1) a careful selection of stars which 
are very similar to each other (stellar twins), and 2) 
a strict differential analysis using spectra of extremely 
high quality. When studying samples of stellar twins, 
the systematic uncertainties that plague chemical abun¬ 
dance analysis are so similar among all the sample stars 
that a strict differential analysis (e.g., measurement of 
relative abundances on a line-by-line basis) essentially 
cancel out. This leaves the observational noise as the 
main source of error. Thus, error bars can be made very 
small simply by acquiring high-resolution spectra of very 
high signal-to-noise ratio {S/N). Instead of the typical 
S/N ^ 100 (or somewhat higher when given per resolu¬ 
tion element) used in standard chemical analyses, all our 
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high-precision abundance studies have employed spectra 
of S/N > 300 per pixel (and typically higher than 400) 
and resolution R = A/AA > 60 000. 

Thanks to the high precision of our derived abun¬ 
dances, we have uncovered interesting abundance anoma- 
lies among stars know n to host planets of different kinds. 


Melendez et al. (2009) first showed that the solar photo- 
sphere is slightly deficient in refractory elements relative 
to volati les when compared to a sam ple of solar twin stars 
(see also Ramfrez et al.|2009 20101. They attributed this 
deficiency to the formation of rocky material in the solar 
system. Essentially, they claim that the missing mass of 
refractories in the Sun was locked-up in the terrestrial 
planets, asteroids, and other solar system rocks at the 
time of star and planet forma tion. 

Later, Ramfrez et al. (2011) showed that the secondary 
star of the 16 Cygni binary system, which is known to 
host a gas giant planet, is slightly metal-poor relative to 
the primary, which has not yet shown evidence of host¬ 
ing planets. Under the assumption that stars in multiple 
systems must have had the same initial chemical compo¬ 
sition, this observation was interpreted as a signature of 
the formation of the gas giant planet, which took both 
refractory and volatile elements from its host star during 
the star/planet formation stage. 

Since the bulk metallicity of gas giant planets is ex¬ 
pected to be higher than that of their parent starsjjthe 
end result of their formation is a decrease of the m etallic- 
ity of the host star. Indeed, Ramfrez et al. (2011) found 
that all elements are deficient by about the same amount 
in 16 Cygni B rela tive to 16 Cygni A. More recently, |Tucci| 


Maia et al. (2014) have shown that there is in fact a very 
small dehciency of refractory elements on top of that con- 


® This is true in our solar system. Using the standard bulk abun¬ 
dance notation (X: h ydrogen, Y: helium. Z: metals), (Z/X) = 
0.04 — 0.12 f or Jupiter iRamirez e t al.|2011|i while {ZfX)Q = 0.018 
for the Sun ijAspiund et aI.|2UUii[l. ' ' 
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slant metallicity offset, which could be the signature of 
the gas giant planet’s rocky core. 

The prospect of detecting and characterizing signa¬ 
tures of planet formation using stellar chemical abun¬ 
dances is tantalizing. Eew other observations allow us to 
probe into these processes, and they tend to be costly, in 
addition to being limited to a handful of very bright and 
nearby young stars. 

In this paper, we investigate the chemical composi¬ 
tion of a binary system known to host planets, XO-2, 
in order to provide further constrains to our hypothesis 
that planet formation imprints signatures in the chem¬ 
ical composition of the host stars. As in our 16 Cygni 
experiment, we begin with the assumption that the two 
stars in this binary system formed from the same gas 
cloud, which had a homogeneous chemical composition. 
Any differences seen in the photospheres of the stars to¬ 
day can then be attributed to processes that occurred 
either during the planet formation stage or at any other 
time during the main-sequence lives of the stars. As far 
as we know, the only previous studies dealing with de¬ 
tailed multi-element chemical abundance ana lyses of this 
system are those by Teske et al. (2013 2015). 


2. THE XO-2 STARS AND THEIR PLANETARY SYSTEMS 

XO-2 (TYC 3413-5-1) is a twin-star binary system 
where planets have been detected around both compo¬ 
nents. The stars are both late G-type or early K-type 
dwarfs of super-solar metallicity. They are separated by 
~ 30arcsec, which corresponds to a projected distance 
of about 4 500AU. The difference in effective tempera¬ 
ture between them is around 60 K, making this system 
comparable to 16 Cygni in terms of the striking similarity 
between the two stars. This gives us confidence that a 
reliable high-precision relative chemical abundance anal- 
y sis can be app l ied to them. 


Burke et al. 


(2007) first announced the discovery of 
a transiting planet around XO-2N (the “North” com¬ 
ponent of the binary). Further follow-up of the system 
has revealed the presenc e of two planets arou nd XO-2S 
(the “South” component; Desidera et al.|2()T4 ) as well as 
evidence of long-t erm variability m the radial velocities 
(RVs) of XO-2N (Damasso et al. 2015, hereafter D15). 
The latter could be interpreted as an additional planet 
around XO-2N, but the detailed analysis of stellar activ¬ 
ity indices by D15 suggests that this variation is likely 
due to the magnetic cycle of XO-2N instead. According 
to these authors, who have performed the most recent 
comprehensive analysis of the XO-2 system and its plan¬ 
ets, the mass of the transiting planet orbiting XO-2N is 
about 0.6 Mj, where Mj is the mass of Jupiter, while 
the masses of the planets orbiting XO-2S, which do not 
transit, are > 0.26Mj and > 1.4Mj. If the long-term RV 
variations detected in XO-2N are due to another planet, 
that object would have a mass > 2.4 Mj. 

Clearly, the planet population around XO-2 is more 
complex than the one around 16 Cygni. In the latter, 
only one planet of mass > 1.5 Mj is known to orbit one 
the stars, namely the secondary 16 Cygni B. The primary 
has not yet shown any eviden ce of substellar mass com¬ 


panions (Cochran et al. 
tation of the observed c 


1997). This made the interpre- 


lemical abundance anomalies in 


16 Cyg ni by Ramirez et al. (2011) and Tucci Maia et al. 
(2014) relatively straightforward. XO-2S has two known 


planets with a total mass > 1.66 Mj. This value could 
be higher or lower than the total mass of planets around 
XO-2N depending on the real nature of the long-term RV 
variability of the latter. And of course we should always 
be cautious about the fact that smaller planets may have 
formed around both of these systems, in one or in both 
of their components, but remain undetected. 

3. SPECTROSCOPIC DATA 

Spectroscopic observations of the XO-2 system were 
carried out with the HIRES spectrograph on Keck’s 10 m 
Telescope on 24-25 January 2015. We used the kv408 
filter, which allows a wavelength coverage from about 
4080 to 8300 A. The slit width was set to 0.57 arcsec, 
corresponding to a spectral resolution R = 67 000. Four 
20 minute exposures for each star were taken, leading to 
a total S/N ~ 350 per pixel in the 5000-7500 A region. 
We also took a solar spectrum of very high S/N (> 700 
per pixel) at R = 84 000. Spectra were reduced using 
the Keck-MAKEE standard pipeline which performs bias 
subtraction, flat-fielding, scattered-light subtraction, ID 
spectral extraction and wavelength calibration. 

Continuum-normalization and merging of spectral or¬ 
ders (hereafter referred to as “apertures”) were per¬ 
formed using IRAF’s tasks continuum and scopy. A 
different polynomial order and pixel sampling was em¬ 
ployed for each aperture in order to visually confirm the 
quality of the continuum-tracing in each star. To ensure 
star-to-star consistency, we adopted the same polynomial 
orders and pixel sampling for every aperture of the three 
(previously co-added) spectra available. 

4. ANALYSIS 


4.1. Input data and analysis tools 

Elemental abundances, including those for iron, which 
are also used to constrain the relative atmospheric pa¬ 
rameters, were determined using a curve-of-growth ap¬ 
proach. Specifically, line strength measurements (equiva¬ 
lent widths) were translated into abundances in the stel¬ 
lar photospheres using standard ID-LTE model atmo¬ 
spheres. The latter were i nterpolated linearly with in the 
Kurucz odfnew grid (e.g., Castelli fc Kurucz‘j|2003 ). The 
abfind driver in the 2014 versio n of the code MOOG 
( Sneden|[l973 Sobeck et al.||2011 ) was employed for the 
spectral line calculations. 

Equivalent widths (EWs) of spectral lines were mea¬ 
sured by fitting Gaussian functions to the observed line 
profiles using IRAF’s splot task. The majority of lines 
in our linelist, described in more detail later, are un¬ 
blended. Those which are somewhat blended typically 
correspond to features due to elements for which only a 
handful of spectral lines are available. 

Our chemical abundance results can be fully repro¬ 
duced using the Qoyllur-quipu (q^) Python package]^ 
Tables provided in this paper contain our EW measure¬ 
ments, and the Keck spectra we employed to make those 
measurements are available upon request. 


4.2. Atmospheric parameters from an iron-line-only 
analysis and using a solar spectrum as reference 

Atmospheric parameters Tgg (effective temperature), 
log 5 (surface gravity), [Fe/H] (iron abundance), and Vt 


https: / / github.com/astroChasqui/q2 
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X0-2N : 5307, 4.3, 0.409, 0.93 [iris] 


XO-2S : 5374, 4.33, 0.351, 0.94 [iris] 
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Figure 1. Iron abundances measured in the spectrum of XO-2N 
relative to the solar abundances as a function of the spectral lines’ 
excitation potential (top panel), reduced equivalent width (middle 
panel), and wavelength (bottom panel). Blue crosses (green circles) 
correspond to Fe I (Fe II) lines. In the top and middle panels, solid 
lines are linear fits to the Fe I data. In the lower panel, the solid 
line is at the average [Fe/H] value. 


(microturbulent velocity) can be measured using only the 
observed spectra. By measuring line-by-line differential 
iron abundances from Fe I and Fe II lines and minimiz¬ 
ing the correlations between relative abundance and ex¬ 
citation potential and reduced equivalent width of the 
lines {REW = \ogEW/\), which can be achieved by it¬ 
eratively fine-tuning the atmospheric parameters them¬ 
selves, it is possible to find a set of values that satisfy the 
so-called excitation and ionization balance conditions. 


This pr ocedure is described in detail in Ramirez et al. 
(2014b), where an explanation of the formal error deter- 
mmation can also be found. In particular, we note that 
the error in the iron abundance inferred is not just the 
1 cr line-to-line abundance scatter obtained with the final 
solution for the parameters, but the result of a rigorous 
error propagation analysis which is provided fo r exam¬ 
ple in Epste in et al. (2010 their Section 3.2) and Bensby 
et al. ([2014 their Appenaix B). 


We employed this method to determine stellar param¬ 
eters for the XO-2 stars using our solar spectrum as ref¬ 
erence, mainly in order to minimize the impact of the un¬ 
certainties in the transition probabilities adopted. As ini¬ 
tial parameters we employed those obtained in the “Dif¬ 
ferential analysis” of D15 (their Tabled, for convenience 


Figure 2. As in Figure[^for XO-2S. 


copied here in the first block of our Tablef^. For the Sun, 
we adopted Teg = 5777 K, logg = 4.44, ^e/H] = 0, and 
Vt = l.lkms“^, as in D15. Our final results, with the 
converged parameters, are shown in Figuresandand 
the derived parameters with their errors are listed in the 
second block of Table[l] The method used to derive these 
parameters ensures that the slopes of [Fe/H] versus exci¬ 
tation potential and REW are zero within error, but it 
does not require a zero slope for the [Fe/H] versus wave¬ 
length relation. Indeed, the latter slopes for Figures 
andj^are not negligible; they are about -1-0.02± 0.01 dex 
per 1000 A. This could be related to systematic effects 
on the continuum fluxes predicted by model atmospheres 
or due to difficulties in measuring equivalent widths in 
the spectra of very metal-rich stars relative to the Sun. 
Indeed, this wavelength trend is closer to zero when the 
XO-2 stars are dire ctly compared to each other, as will 
be shown in Section 

Figures [T and [^illustrate some of the most important 
features ofour analysis. Note the large number of Fe I 
and in particular Fe II features employed. Note also that 
the excitation potential coverage of our Fe I linelist is 
such that an important number of low excitation poten¬ 
tial (EP) Fe I features are available. These lines are very 
important for constraining the effective temperature of 
the star within the iron-line-only analysis because the 
EP trend is highly sensitive to Teff. Other works tend 
to have few of those lines because they are difficult to 
measure (many low EP Fe I lines require very careful de- 
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Table 1 

Atmospheric Parameters 


Star 


'Teff 

(K) 


logs 

[£Si 


[Fe/H] 
(dex) 


vt 

(kms“^) 


XO-2N 

XO-2S 

A(N-S) 


5290±18 
5325 ±37 4.42 ±0.10 

-35 ±8 ±0.01 ±0.02 


Damasso et al.|(|2015 I 

4.43 ±0.10 0.37 ± 0.07 


0.32 ±0.08 
±0.054 ±0.013 


0.86 ± 0.06 
0.93 ± 0.03 
-0.07 ± 0.07 


Iron-line-only with solar reference (Section 

XO-2N 5307 ± 19 4.30 ± 0.05 0.41 ± 0.02 

XO-2S 5374 ± 16 4.33 ± 0.04 0.35 ± 0.02 

A(N-S) -67 ± 25 -0.03 ±0.06 ±0.058 ± 0.023 


4.21 


0.93 ± 0.05 
0.94 ± 0.04 
-0.01 ± 0.06 


Strict N— S comparison, IRFM Tgg and “high” logg (Section 4.51 
XO-2N 5440 4.43 0.40 0.94 

XO-2S 5500 ±5 4.45 ±0.02 0.347 ±0.004 0.94 ± 0.02 

A(N-S) -60 ±5 -0.02 ±0.02 ±0.053 ± 0.004 0.00 ± 0.02 


blending). Instead of rejecting those lines, we make an 
effort to measure as many of them as precisely as possi¬ 
ble. 

The parameters we derived as described in this section 
are in reasonably good agreement with those by D15. 
Our Teff values are slightly warmer, more so for XO-2S 
than XO-2N, which implies a larger N-S T^s difference 
in our case (—67 instead of —35K). Also, while D15 sug¬ 
gest a logg ~ 4.42 for both stars, we find lower values 
(logg ~ 4.32). This discrepancy is in principle consis¬ 
tent within the error bars for each individual star, but 
since we detect about the same offset for both stars, it 
is unlikely that this difference is just a product of ob¬ 
servational noise. Although we do not require extremely 
accurate logg values in order to measure very precise 
relative abundances, they are important for determining 
the evolu tionary state and age of the stars. Thus, in Sec¬ 
tion [4^ we continue our discussion of the logg values. 

Eigures[l] and show that the iron abundance of XO- 
2N is slightly higher than that of XO-2S. The iron abun¬ 
dance difference (N—S) we find with this procedure is 
±0.06 ± 0.02, in good agreement with D15: ±0.05 ± 0.01 
(±0.05 ± 0.11 when using a solar spectrum as reference). 


4.3. Photometric effective temperatures 

Effective temperatures of stars can also be measured 
using photometric methods. Since they rely on broad¬ 
band features rather than spectral lines, they are ex¬ 
pected to be less sensitive to model uncertainties than 
iron-line-only analyses. One of the most powerful pho¬ 
tometric effective temperature determination techniques 
is the so-called InfraRed Flux Method (IRFM ; e.g. , 
Blackwell & Shallis 1977[ Ramirez fc Melendez 2005 


Casagrande et al .||2066 2(jl0p. 

Using the photometric data provided in Table 1 of the 
D15 paper, and adopting their E{B — V) = 0.019 value, 
we computed IRFM effective temperatures for the XO-2 
stars assuming [Fe/H] = +0.4 and logg = 4.5 for both 
stars. The IRFM is very weakly dependent on [Fe/H] and 
logg so our use of these approximate values is accept- 
able. The IRFM i mplem entation we adopted is that by 
Casagrande et al.|(|2010|. In this scheme, three tempera- 
ture values are determined from the J, H, and Kg bands, 
respectively. The standard deviation of these three val¬ 
ues gives us a first estimate of the error. To account 
for conservative uncertainties of about 0.1 dex in [Fe/H] 
and 0.005 mag in E{B — V), we added in quadrature 


Table 2 

Photometric Effective Temperatures 


Star 

Direct IRFM 
(K) 

Color calibrations 
(K) 

Weighted average 
(K) 

XO-2N 

5472 ± 48 

5414± 56 

5447 ± 36 

XO-2S 

5530 ± 42 

5474 ± 50 

5507 ±32 

A(N-S) 

-58 ± 64 

-60 ± 75 

-60 ± 48 

A(N-S) 

— 

-57± iFl 

— 


^Using only the differential photometry by [Damasso et ah](|2015]| 


errors of 10 K and 30 K, respectively, to the IRFM 
error. The Teff values obtained using the IRFM directly, 
and their errors, are listed in the second column of Ta- 
bl^( “Direct IRFM”). 

Tne stars’ colors can also be employed to measure effec¬ 
tive tempe ratures by means of Teff-c olor calibrations such 
as those by Casagrande et al. (2010). Since these calibra¬ 
tions are based on IRh'M 'j/ff values determined in the 
same manner we calculated the IRFM temperatures of 
the XO-2 stars above, using Teff-color calibrations allows 
us to reduce the impact of errors in the photometric data 
on the final, averaged Teff values, while keeping them con¬ 
sistent with the same IRFM effective temperature scale. 
The latter has been shown to be in good agreement with 
Teff scales based on direct measurements of angular di¬ 
ameters and bolometric fluxes of a selected number of 
nearby dwarf and sub giant stars (Casagrande et al. ]mol 
Boyajian et al.| 2013), so it is desirable to remain consis¬ 


tent with that 'Jeff scale. We calculated {B — V), (V—J), 
{V — El), and {V — Kg) colors for the XO-2 stars using 
the photometric data from Table 1 in D15, and adopt- 
ing the extinction ratios, k = E{co\ot) /E{B — V), from 


Ramirez & Melendez (2005). The average Teff values ob- 
tained using IRh'M Jeff-color calibrations are listed in 
the third column of Table [^ In this table we also pro¬ 
vide the weighted averages ot direct-IRFM and Teff-color 
values, which we adopt as our final photometric effective 
temperatures (fourth column of Table 
The photometric temperatures for the XO-2 stars are 
significantly warmer (about +140 K) than those obtained 
in the iron-line-only analysis. The difference is non- 
negligible, not within the 1-cr error bars, and it has the 
same sign for both stars, suggesting that it is due to 
systematic errors. Although the photometric Teff values 
that we have determined are warmer than those we mea¬ 
sured with our iron-line-only analysis as well as those 
from the D15 work, we should note that they are not 





















































XO-2 ELEMENTAL ABUNDANCES 


5 


too far off from the v alues given in Torres et al. (2012), 


Mortier et al.| (|2013|) , and |Teske et al.| (j2U15(), for ex¬ 
ample. i’he tormer used spectrum matcning to a grid 
of synthetic spectra while the latter two employed iron- 
line-only techniques similar to ours. 

According to the photometry, the Tes difference be¬ 
tween the XO-2 stars (N-S) is —60 ± 48 K. Although in 
principle consistent with the difference reported by D15 
(—35K) within the error bars, the mean value we obtain 
from our iron-line only analysis is closer (—67K). More¬ 
over, using the very precise differential photometry pro¬ 
vided by D15 (given in the second block of their Table 1) 
we can also infer a photometric Teg difference which is 
nearly insensitive to the exact value of the stars’ colors. 
The differential {B — V), {V — R), and (U —/) colors have 
such small errors (0.002-0.004 mag) that they translate 
into photometric T^ff difference errors of less than about 
10 K. Averaging the three differential Tes values obtained 
in this manner we determined a much more precise pho¬ 
tometric N-S Teff difference of —57 ± 7K. 

4.4. Surface gravity and evolutionary state 
As pointed out in Section |4.2[ the logg value that 


we determined from our iron-line-only analysis is about 
0.1 dex lower than that obtained by D15. To place this 
important atmospheric parameter in context, in Figure!^ 
we show theoretical isochrones of 1, 5, and 9Gyr (and 
[Fe/H] = -1-0.4) on the Tes-logg plane along with the 
location of the XO-2 stars depending on the method of 
atmospheric parameter determination used. Our “low” 
log g values would make both stars fall outside of the al¬ 
lowed region for Galactic disk stars, i.e., older than about 
8 Gyr (red open symbols). Similarly, the parameters by 
D15 (blue half-filled symbols) would make the system a 
very old one (> 9Gyr). 

If the effective temperatures of the XO-2 stars are 
closer to 5500 K, as implied by the photometry, the range 
of allowed log 5 values expands to ~ 4.3 — 4.5. Nev¬ 
ertheless, when we fix the Teff of the XO-2 stars to 
our photometrically-determined values, therefore basi¬ 
cally just forcing the Fe I/Fe II ionization equilibrium, 
we find that we require logg ~ 4.6 for both stars, which 
is clearly inconsistent with the predictions of stellar evo¬ 
lution theory. The low levels of chromospheric activ¬ 
ity detected in these stars and their long rotation peri¬ 
ods suggest a moderately old age. However, the stars 
are still members of the Galactic thin disk, as argued 
by D15. Therefore, they must be younger than about 
8 Gyr. Indeed, using the rotation period of XO-2N ac¬ 
curately measured by D15, which is about 42 days, and 
its {B — V) color, the age ( pf the system according to the 
gyrochronology formula of Barnes (2007) is ~ 6 Gyr. For 
this age (and Teg ~ 5500K), the preferred logg values 
for the XO-2 stars are ~ 4.45. By fixing the and 
logg values of XO-2S as suggested above, we deri ved pa¬ 
rameters for XO-2N as described in Section 14.51 These 
parameters are shown in Figure!^ as green solid symbols. 

The parameters from D15 result in an excellent agree¬ 
ment with the isochrones, i.e., both stars fall almost per¬ 
fectly on the 9 Gyr Yonsei-Yale isochrone, as expected 
for a coeval system of two stars. On the other hand, our 
preferred mean parameters, the green solid symbols, do 
not fall exactly on the same isochrone. However, con¬ 
sidering the uncertainties in absolute stellar parameters 



Figure 3. Isochrones of 1, 5, and 9 Gyr as predicted by Yonsei- 
Yale (dashed lines) and Padova (solid lines) stellar evolution models 
in the Teff-logp- plane. Squares and circles represent the location 
of the XO-2S and XO-2N stars, respectively, for three different 
choices of stellar parameter determination: iron-line-only (red open 
symbols), photometric Teff and lo g q = 4.45 fixed for XC )-2S (green 
filled symbols), and the values by|Damasso et al.|(|2015|l (blue half- 
filled symbols). Representative error bars tor ttie absolute values 
of Teff and log^ are shown for the data point at Teff = 5500 K, 
logg = 4.45. 

(about 30 K in Teff from the IRFM and about 0.04 dex in 
logg from the iron-line-only analysis relative to the Sun), 
our preferred values are reasonably consistent with the 
expectation based on predictions of standard stellar evo¬ 
lution theory. 

4.5. Parameters from a strict N—S comparison 

Spectroscopic equilibrium, or excitation/ionization 
balance of iron lines, is far from being satisfied when em¬ 
ploying photometric Teft values and a “high” logg ~ 4.45 
for the XO-2 stars. The trends with excitation potential 
and reduced equivalent width are highly significant, as 
is the difference between iron abundances inferred from 
Fe I and Fe II lines. These trends likely reflect the lim¬ 
itations of traditional ID-LTE analyses. Future studies 
with tailored 3D models and accounting for non-LTE ef¬ 
fects are of great interest, but beyond the scope of this 
paper. From a more practical point of view, one can try 
to avoid these issues by not comparing the XO-2 stars 
to the Sun anymore and performing a strict differential 
analysis of XO-2N relative to XO-2S instead, after adopt¬ 
ing a set of parameters for XO-2S. 

The stellar parameters that we adopt for XO-2S are 
not entirely arbitrary. Photometric effective tempera¬ 
tures and surface gravities consistent with stellar evolu¬ 
tion predictions are more reliable in an absolute sense 
than Teff, logg combinations that result from purely 
spectroscopic analyses of iron lines. Thus, hereafter we 
fix Teff = 5500K, logg = 4.45, [Fe/H] = -tO.35, and 
Vt = 0.94 for XO-2S, and measure only relative parame¬ 
ters for XO-2N. The solar spectrum is no longer used to 
calculate stellar parameters. The technique employed to 
derive the relative parameters is the same one used be¬ 
fore with the only exception that XO-2S is the reference 
star instead of the Sun. 

The final result of our strict line-by-line differential 
stellar parameter determination is shown in Figure 
The derived parameters for XO-2N and their formal er¬ 
rors are listed in the third block of Table [H Note that 
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X0-2N : 5440, 4.41, 0.4, 0.94 [XO-2S] 
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Figure 4. As in Figure for the line-by-line analysis of XO-2N 
relative to XO-2S. Note tnat the scale of the vertical axis here is 
about four times smaller than that in Figures]^ andand that the 
line-to-line scatter here is about six times smmler because the two 
stars being compared are very similar to each other, minimizing 
systematic errors. 

the latter are considerably smaller than those obtained 
before. This is due to the fact that we are now comparing 
two stars which are twins of each other, and do not use a 
solar spectrum as an intermediate stepQ In Sectionj^ we 
advocated using the Sun as a comparison. This approach 
is of great value when analyzing solar twins. Comparing 
abundances of XO-2N and XO-2S with each other rather 
than the Sun significantly red uces systematic errors. In¬ 
deed, as noted in Section [T^ a non-negligible slope for 
the [Fe/H] versus wavelength data is observed when com¬ 
paring these stars to the Sun. The corresponding slope 
for Figure 0 (A[Ee/H] versus wavelength) is about four 
times smaller: -1-0.005 ± 0.002 dex per 1000 A, yet still 
not zero within the very small errors in this twin-star 
comparison. 

Our strict N-S comparison, in which we adopted a 
photometric (IRFM) effective temperature and a “high” 

^ Since the EWs for the iron lines were measured “manually” by 
one of us, instead of using an automated code, one could argue that 
our results are subjective. To address this potential problem, an¬ 
other set of iron line EtFs were measured by a different person and 
the differential parameters were determined using this independent 
set of EWs. In that case we found ATgff = —65 ± 8K, Alogg = 
-0.03 ± 0.03, A[Fe/H] = -1-0.063 ± 0.010, and Avt = -0.03 ± 0.03, 
in very good agreement with the first case, albeit with somewhat 
larger error bars. This shows that with very high-quality spectra 
the subjectivity of manual EW measurements is not so important. 
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surface gravity for the reference star XO-2S, reveals a 
very precise effective temperature difference (N-S) of 
—60 ± 5 K, in excellent agreement with the very pre¬ 
cise value obtained from the differential photometry 
(—57 ± 7K). Both are significantly higher than the dif¬ 
ference measured by D15 (—35 ± 8 K). 

Interestingly, despite the differences in stellar parame¬ 
ters adopted and the Tes difference, the offset in [Fe/H] 
between XO-2N and XO-2S is almost exactly the same: 
-1-0.054 ± 0.013 according to D15 versus our -1-0.053 ± 
0.004. In fact, any of the other combinations of stellar 
parameters that we have found so far always result in 
an iron-enhanced XO-2N relative to XO-2S by at least 
-1-0.05 dex. There is little doubt that the iron abundances 
of these two stars are different. 

According to our strict N-S comparison, A log g = 
—0.02±0.02 (N-S). This means that the cooler N compo¬ 
nent has a slightly lower log g, which implies an isochrone 
age about 1 Gyr older than the S component, in contra¬ 
diction with the assumption of coeval age for stars in 
binary systems. Note, however, that the difference is 
marginally consistent with zer o, and in fact a deta iled 
isochrone age calculation (as in Ramfrez et al. |2014bD re¬ 
sults in error bars of about 1 Gyr tor each star. INeverthe- 
less, it is interesting to note that adopting a log <7 = 4.45 
for XO-2N results in an Fe II minus Fe I iron abundance 
difference of -bO.Ol dex, as opposed to exactly zero. This 
slightly higher Fe II abundance could be related to the 
Fe II minus Fe I discrep ancies seen in the Hyades cluster 
by Yong et al. (2004), who employed log 5 values de- 
termined from isochrones and a photometric scale. 
Thus, it is possible that the “peculiar” log g value of XO- 
2N is again a product of model deficiencies. However, as 
we discuss later in this paper, this issue does not affect 
our relative abundance ratios in any significant way. 

4.6. Elemental abundance determination 

Elemental abundances of 22 elements other than iron 
were measured for the XO-2 stars. For carbon, we used 
GI and CH lines, while for Sc, Ti, and Gr, spectral lines 
due to the neutral and singly-ionized species were em¬ 
ployed. The li nelist and atomic para meters adopted are 
the same as in Ramfrez et al. (2014b), except for the GH 
and Rb lines, t^r these species we adopt ed the line data 
compi led by Asplund et al. (2005) and Grevesse et al. 
(2015), respectively. Hypernne structure was taken into 


account for V, Mn, Go, Gu, Rb, Y, and Ba. The linelist 
and atomic parameters used in this work are listed in 
Table along with the EWs measured in the XO-2 stars 
and our solar spectrum. The relative abundances ob¬ 
tained from each species, and their errors, are given in 
Table The oxygen abundances inferred from the O i 
777 nm triplet were corrected for differential non-LTE ef- 
fects using the grid and code by [Ramfrez et al.| (|2007|). 

For each species, we calculated average values of 
A[X/H](N-S) and the standard deviation of the line- 
to-line scatter (cr), from which a standard error was 
computed as dimes = cr/Vn — 1 , where n is the num¬ 
ber of spectral lines employed. In cases where only 
one spectral line is available we conservatively adopted 
diines = 0.025 dex, which corresponds to the largest error 
obtained for species with more than three lines available. 
On the other hand, by propagating the formal errors in 
the relative atmospheric parameters of XO-2N, as given 
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Table 3 

Line List 


Wavelength 

(A) 

Species 

EP 

(eV) 

logs/ 

XO-2N 

(mA) 

XO-2S 

(mA) 

Sun 

(mA) 

5052.167 

6.0 

7.685 

-1.304 

32.0 

34.2 

32.8 

6587.6099 

6.0 

8.537 

-1.021 

10.8 

12.1 

13.8 

7111.469 

6.0 

8.64 

-1.074 

17.7 

19.1 

10.7 

7116.96 

6.0 

8.65 

-0.91 

11.9 

12.7 

13.2 

4218.723 

106.0 

0.413 

-1.008 

92.5 

93.9 

79.5 

4253.003 

106.0 

0.523 

-1.523 

43.3 

44.3 

33.3 

4253.209 

106.0 

0.523 

-1.486 

42.2 

43.8 

33.8 

7771.9438 

8.0 

9.146 

0.352 

44.2 

48.3 

68.7 

7775.3901 

8.0 

9.146 

0.002 

38.9 

42.0 

48.0 

4751.8218 

11.0 

2.104 

-2.078 

39.1 

34.9 

18.4 

4982.82 

11.0 

2.1 

-1.0 

103.5 

101.9 

81.9 

5148.838 

11.0 

2.102 

-2.044 

40.9 

36.3 

10.0 


Table 4 

Relative Abundances 


Species 

7c 

(K) 

A[X/H] 

(dex) 

error 

(dex) 

Cl 

40 

0.009 

0.009 

CH 

40 

0.015 

0.006 

Oi 

180 

0.013 

0.008 

Nai 

958 

0.035 

0.009 

Mg I 

1336 

0.074 

0.012 

All 

1653 

0.087 

0.018 

Si I 

1310 

0.070 

0.006 

Si 

664 

0.024 

0.008 

Ki 

1006 

0.008 

0.027 

Cai 

1517 

0.071 

0.009 

Sc I 

1659 

0.102 

0.024 

Sen 

1659 

0.076 

0.010 

Til 

1582 

0.085 

0.009 

Ti II 

1582 

0.082 

0.010 

Vi 

1429 

0.078 

0.009 

Cri 

1296 

0.044 

0.008 

Cr II 

1296 

0.067 

0.011 

Mnl 

1158 

0.035 

0.008 

Fel 

1334 

0.054 

0.005 

Fell 

1334 

0.058 

0.007 

Co I 

1352 

0.058 

0.006 

Nil 

1353 

0.061 

0.006 

Cui 

1037 

0.055 

0.011 

Zn I 

726 

0.020 

0.006 

Rbi 

800 

0.019 

0.017 

Yii 

1659 

0.082 

0.019 

Zr II 

1741 

0.116 

0.026 

Ball 

1455 

0.074 

0.017 


trend perfectly (Tc = 1334 K for iron)j^ 

A double linear fit to the A[X/H] versus Tq data is 
also shown in Fig ure We used the IDL MPFIT routine 
(Markwardt 2009) to nnd a robust value for the Tc of the 
break. This routine uses the Levenberg-Marquardt tech¬ 
nique for least-squares minimization and requires an ini¬ 
tial guess. We tested a large range of initial guesses and 
found the break to be at Tq = 841 ± 137 K. With a break 
at Tc = 841K the connection between the two linear fits 
is smooth. Each linear fit was calculated using the inverse 
square of the A[X/H] errors as weights. The minimum 
and maximum A[X/H] values of the fit are 0.014 and 
0.090 dex, respectively. The element-to-element scatter 
around the fit is 0.011 dex, which is very similar to the 
average value of our A[X/H] errors (0.012 dex). This sug¬ 
gests that the linear fits shown in Figure [^constitute an 
excellent representation of the abundanceaifferences be¬ 
tween the XO-2N and XO-2S stars and that the element- 
to-element scatter is due to observational errors alone. 
Furthermore, this implies that the A[X/H] difference for 
the most volatile elements such as C and O is slightly 
larger than zero. In fact the volatile element difference 
can be quantified as -1-0.015 ± 0.011 dex. Since the most 
refractory elements are over-abundant in XO-2N relative 
to XO-2S by about 0.090 dex, the amplitude of the Tc 
trend shown in Figure is about 0.075 dex. 


4.7. Dependence on adopted stellar parameters 

Here we discuss how the result in Figure depends 
on alternative choices of stellar parameters. Through¬ 
out this exercise, we have used our EW measurements. 
Thus, it is not our intention to replicate the results of 
previous studies, whenever applicable, but rather to test 
whether the very strong A[X/H] versus Tc correlation 
that we have found is only due to our particular choice 
of stellar parameters or their N-S difference. 

We begin by recalculating the abundances using our 
iron-line-only parameters from Section |4.2[ which em¬ 
ploys the Sun as the reference star. These parameters 
are given in the middle section of Table The most 
important differences with our preferred parameters are 
the lower Teff and logg values, but the N-S differences 
are essentially the same. Panel a) in Figure]^ shows the 
XO-2N minus XO-2S abundance difference m this case. 
The Tc trend is still very clear and there appears to be 
only a minor offset of about —0.01 dex relative to our pre¬ 
ferred abundance set. In this case the volatile element 
abundances are slightly more similar between XO-2N and 
XO-2S, but the latter still seems to be slightly volatile 


in the last row of Table we computed an uncertainty 
in the average A[X/H] values, dpars, which varies from 
one species to another in the 0.003-0.011 dex range. We 
added (funes and ^pars in quadrature to obtain the final 
error bars for our A[X/H](N-S) relative abundances. 

Figure^ shows the elemental abundance difference be¬ 
tween XO-2N and XO-2S as a functio n of condensa- 
tion temperature, the latter as given by Loddersj (2003 


specifically their 50 % condensation sequence tor a solar- 
composition mixture). The correlation between A[X/H] 
and condensation temperature is very strong while the 
~ -1-0.05 dex offset in the iron abundance previously re¬ 
ported by D15 and confirmed in this work fits the Tc 


® Most, if not all studies that have looked for correlations be- 
tween abundance anomalies and condensation temperature have 
employed the solar Tq sequence as reference (i.e., as the x-axis 
in figures like our Figure and for statistics associated with it. 
However, the Tc values depend on the composition of the gas, 
and one would expect the differences with the solar case to be 
particularly important for metal-rich systems like XO-2. Detailed 
chemical equilibrium calculations to dete rmine the Tn sequ ences 
on a star-by-star basis, such as those by |Bond et aL| ||2010|l, are 
beyond the s cope of this paper. We note, nowever, that tne Tc 
sequences by [Bond et al.| for stars somewhat similar to those in 
the XO-2 system result in noisier Tc trends. Although the overall 
correlation (the “up ward” trend seen in Figure]^ persists, the Tc 
values computed by [Bond et al.| for C and O can be significantly 
higher at high metallicity (up to about 1 000 K) while those of some 
refractories like Si and Ti could be higher by about 300 K. These 
effects should be investigated in future works in this area. 























RAMIREZ ET AL. 


XO-2 (N-S) 



Figure 5. Chemical composition difference between XO-2N and XO-2S as a function of the elements’ condensation temperature. The 
dashed line corresponds to identical chemical composition. The solid line is a double linear fit to the data, broken at Tq = 841 K. Red 
vertical lines connect two species of the same chemical element (e.g., CH and Ci, Til and Tin, etc.). 


poor. Although the values used in this case are about 
140 K cooler than the more reliable IREM t emp eratures, 
and the logg values are too low (see Section 4.4), the fact 
that the N-S relative parameters are basically the same 
leads to a very similar A[X/H] versus Tq correlation. 

Panel b) in Eigure ra sh ows the abundances we derive 
using the parameters from ||Damasso_et al. (2015), which 
are given in the first section of 'i'able[2 Perhaps the most 
important difference with our parameters is the ATeff 
value, which is about twice as large in our case (D15 
suggest ATeff = —35K, but we find ATeff = —60K). 
Also, the Teff values themselves are offset by about 160 K, 
although this might be of secondary importance consid¬ 
ering the test discussed in the previous paragraph. 

Using the parameters from D15, the A[X/H] versus Tq 
correlation appears steeper, but it is also more noisy than 
the one we obtained using our parameters. We note also 
that in this case certain elements with more than one 
species available for analysis show abundance differences 
which are less consistent than with our parameters. Por 
example the carbon abundances obtained from C i and 
CH lines show a difference of 0.04 dex whereas with our 
parameters the difference is only 0.006 dex. Also, the Ti 
abundance from Ti i lines is about 0.03 dex higher than 
that from Ti ii lines in this case, but we find a difference 
of 0.003 dex with our parameters. Based on these re¬ 
sults, we argue that our differential parameters are more 
accurate than those given by D15. In any case, the Tq 
trend of the elemental abundance difference is not at all 
removed when using their parameters. In fact, as noted 
above, the Tc trend seems to be steeper, with the volatile 


abundances being nearly the same between XO-2N and 
XO-2S, or even slightly lower, but only if the carbon 
abundance from the CH lines is ignored. 

In panel c) of Figure ^ we show the relative abun¬ 
dances that we derive wRen we adopt MARCS model 
atmospheres instead of those from Kurucz’s odfnew grid. 
The differences are all within 0.01 dex. Therefore, adopt¬ 
ing a different model atmosphere grid has no impact on 
the A[X/H] versus Tq correlation found . As already 


shown in some of our previou s papers (e.g., Ramfrez et al. 
2011[ Melendez et al. 2012), the choice of model atmo¬ 


sphere grid is irrelevant when dealing with high-precision 
(~ 0.01 dex) relative abundances of stellar twins. 

Based on the discussion above, it appears to be the 
case that a very precise knowledge of the temperature 
difference between the XO-2 stars is critical to trace the 
A[X/H] versus Tc correlation. This would be expected if 
the average excitation potential of the lines employed for 
each chemical element correlates with the condensation 
temperature. Figure [^shows that this is in fact the case. 
With few exceptions (CH, Rb, and K), the average EP of 
the lines used for a given element is correlated with the 
elements’ Tq. Thus, one might wonder what ATeff value 
would completely blur the A[X/H] versus Tq correlation 
and/or flip the sign of the slope. 

Figure shows the results obtained for the XO-2N 
minus X(>2S abundance difference for various arbitrary 
choices of ATeff (N-S). Only the effective temperature dif¬ 
ference is changed; all input data as well as the other pa¬ 
rameters, including Tes = 5500 K for XO-2S, were kept 
constant. We varied ATeg from -1-60 K to — 180K, but 
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Figure 6. As in Figure!^ but for different choices of stellar pa¬ 
rameters for XO-2N and XO-2S. The solid lines correspond to the 
fit to our data in Figure]^ 



Figure 7. Average excitation potential of the spectral lines used 
in this work as a function of the elements’ condensation temper¬ 
atures. The three species which depart the most from the mean 
trend (CH, Rb, and K) are labeled. 
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Figure 8. As in Figure but for different values of the XO- 
2N minus XO-2S temperature difference, ATeff, keeping all other 
parameters (from the third block of Table [y constant. Note that 
panel c) is a replica of Figure The solid nnes correspond to the 
fit to panel c)’s data. The three species that depart the most from 
the average EP versus Tq correlation (CH, Rb, and K), shown in 
Figure are labeled. 
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note that positive ATeff values are unrealistic given the 
very precise differential photometry of these objects and 
the fact that they are co-eval main-sequence stars. XO- 
2S is definitely brighter and bluer, therefore warmer than 
XO-2N, which implies ATeff < 0- 

The element-to-element scatter in A[X/H] versus Tq 
increases for both lower and higher values of ATeff rela¬ 
tive to our preferred value (—60 K), which corresponds to 
panel c) in Figure Moreover, the agreement in abun¬ 
dance difference for elements with two species available 
is excellent when using ATgff = —60K, but it deterio¬ 
rates for other values (see the filled circles connected by 
red vertical lines in Figure^. This is particularly the 
case of element C, for whichmigh excitation Ci and low 
excitation CH lines are available. While the C abun¬ 
dance difference from C i lines becomes less negative and 
eventually more positive for higher values of ATgff, the 
C abundance from CH lines shown the opposite behav¬ 
ior. The C abundance difference from both species is 
virtually the same for ATgg = —60K. 

Inspection of Figure shows that the A[X/H] ver¬ 
sus Tc slope remains positive, albeit with progressively 
larger error, up to ATgff = -1-60 K and down to ATgff = 
— I20K, but it changes sign at about ATgg = —I80K, 
particularly if CH lines are not used. Thus, only if the 
absolute value of the temperature difference between XO- 
2N and XO-2S is highly overestimated, by more than 
100 K, would the observed Tc slope change sign. 

Figure could also be used in the XO-2 context to 
explore the possibility that stellar activity affects the de¬ 
rived abundances at high precision. DI5 have noticed 
that in the timespan of their observations XO-2N was on 
average more active than XO-2S. Low excitation lines, 
thought to form preferentially at high atmospheric layers, 
would be more affected by activity-related heating than 
high excitation lines. If this had a significant impact on 
the line strengths, one would expect it to introduce an 
abundance trend with EP (and therefore Tc). 

However, we must note that line formation de pth is 
much more sensitive to line s trength than EP (e.g., Gur- 
tovenko & Sheminova 20151 and that our linelist is ex- 
tremely heterogeneous in this regard. In addition, while 
the different activity levels of the XO-2 stars could be 
explained by different current phases of their activity cy¬ 
cles, one should notice that the observed ranges of ac¬ 
tivity indices in XO-2N and XO-2S overlap significantly 
(Figures 15 and 17 in DI5). Thus, the probability that 
the activity levels of the two stars are very similar at any 
given time of observation is not low. Unfortunately, our 
spectra do not include the Ca II H&K lines. 


4.8. Using parameters from Teske et al. (2015) 


Teske et al. (2015 hereafter TI5) provided an in-depth 
investigation of stellar parameter dependency for the 
chemical abundance difference between the stars in the 
XO-2 system. Four sets of parameters were determined 
by them (they are listed in their Table I), from which 
they calculated abundances of 16 elements using spectra 
of i? ~ 60 000 and S/N ^ 170 — 230 taken with the High 
Dispersion Spectrograph at the 8.2 m Subaru Telescope. 

Figure|^shows the A[X/H] versus Tc relations we find 
when we use our EW measurements, but the stellar pa¬ 
rameters from TI5. We show only the data points that 
correspond to the species used in TI5. Note in particular 


that their C abundances are based on C i lines only. Also, 
while their O abundances were obtained via spectrum 
synthesis of the forbidden 630 nm line, ours were calcu¬ 
lated from EW measurements the O i triplet at 777 nm. 

When using TI5’s “original parameters,” the A[X/H] 
versus Tq relation shows a negative slope, as shown in 
panel a) of Figure Note that in this case two values 
of each Fe and Ti abundances are available, and they 
disagree badly (see the red vertical solid lines). Panel a) 
of Figure [^closely resembles panel e) of Figurej^ Indeed, 
the ATgff^-S) according to TI5’s original parameters is 
—204K, much cooler than our precise value of —60K. In 
any case, we note that our N-S abundance differences for 
this set of parameters are in very good agreement with 
those found by TI5. The general appearance of the black 
asterisks in their Figure I closely resembles that of our 
filled circles in Figure!^). In particular, we note that the 
A[X/H] values for C, U, S, and Zn (all with Tq < 800K), 
are consistently all positive. 

Panels b), c), and d) of Figure show the A[X/H] 
versus Tc correlation that we obtain using TI5’s “al¬ 
ternative parameters” I, 2, and 3, respectively. All of 
them exhibit somewhat steeper slopes compared to our 
mean trend (shown with green solid lines), but also more 
element-to-element scatter. The alternative parameters 

1 and 2 have ATeff = —83K and —72K, values which 
are not very different from ours. For the alternative pa¬ 
rameters set 3, ATeff = +3K. 

While the steeper Tc slope can be explained for the 
alternative parameter set 3 by comparing it to our test 
in panel b) of Figure]^ the steeper slopes and noisier be¬ 
haviors obtained witn the alternative parameters I and 

2 are more likely due to discrepancies in the other stel¬ 
lar parameters. In any case, it is reassuring to find that 
these Tc slopes are positive and that the detailed A[X/H] 
versus Tc relations that we find using our EW measure¬ 
ments are very similar to those obtained by TI5. Indeed, 
the red and blue symbols in their Figure I look remark¬ 
ably similar to our Figure |^s panels b), c), and d). In 
all these cases, both TI5 and ourselves find C, O, S, and 
Zn abundance differences (N-S) below zero (these are 
the four elements with the lowest Tc in Figure]^. At 
the same time, virtually all other species show positive 
differences in both our calculations and TI5’s. In addi¬ 
tion, the average difference for the most refractory ele¬ 
ments available is about -l-O.Idex, which combined with 
the negative differences seen in the volatiles leads to the 
more steep Tq slopes. 

The A[X/H](N-^ versus Tc correlation t hat we have 
found (our Figure™ was already seen in the Teske et al. 


(20151 data, but with more element-to-element scatter. 
We believe that part of this scatter is due to the less pre¬ 
cise parameters within the “Alt Params” sets, but also 
due to the lower quality of their spectra, which trans¬ 
lates to less precise EW values and therefore less accu¬ 
rate relative abundances. On the other hand, it is clear, 
based on our various tests, that the “Original Params” 
set from TI5 does not represent correctly the fundamen¬ 
tal properties of the XO-2 stars and that the abundances 
derived using those parameters are unreliable. Neverthe¬ 
less, their key finding that “Fe, Si, and potentially Ni are 
consistently enhanced in XO-2N” is fully supported by 
our study. 
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Fig ure 9. As in Figure j^ but using input atmospheric parameters 
from|Teske et al.|(|2015|l and restricted to the species used in thid 
work! ihe solid lines correspond to the fit to our data in Figure]^ 

The tests made in this section highlight that our rel¬ 
ative parameters, and very likely their absolute values, 
are more precise than those previously published for the 
XO-2 stars. We attribute this mainly to the higher qual¬ 
ity of our data, but also to our use of more reliable ef¬ 
fective temperature and surface gravity indicators, even 
when they lead to inconsistencies in iron-line-only types 
of analyses. One must realize and accept the fact that 
those inconsistencies are most likely related to our limi¬ 
tations in the modeling of stellar atmospheres and spec¬ 
tral line calculations, and that they should not be used 
as arguments against other techniques of stellar parame¬ 
ter determination, especially when these other techniques 
have proven to be less affected by modeling errors. 


5. PLANET SIGNATURES 
5.1. Metal depletion from planet formation 


Acc ording to Melendez et al. (2009) and Ramirez et al. 
(2011), planet formation results m a small dehcit of met¬ 
als m the host star’s convective envelope relative to its 
initial composition (and also relative to the composition 
of the star’s radiative zone). The formation of a planet 
slightly decreases the overall metallicity of the remain¬ 


ing gas in the proto-stellar nebula. If that gas is then 
accreted onto the host star, its surface metallicity is low¬ 
ered. Specifically, these studies suggest that when giant 
planets form, the overall metallicity is changed, while 
terrestrial planets have a more important effect on the 
abundances of refractory elements compared to volatiles. 

As mentioned in Section!^ a O.GMjup planet has been 
detected by transit and RV aata around XO-2N while two 
planets with minimum masses 0.26 Mjup and 1.4Mjup 
have been inferred from the RV data of XO-2S. In ad¬ 
dition, a long-term RV variability around XO-2N could 
be interpreted as an additional l.SMj^p planet, but its 
strong correlation with stellar activity indicators suggests 
that it is more likely due to a magnetic cycle instead. 

Figure shows that XO-2S is metal-poor relative to 
XO-2N. Simply put, one could say that this result im¬ 
plies that more disk material was used to make planets in 
XO-2S. This would be compatible with the known planet 
population of the system if the long-term RV variability 
of XO-2N is due to magnetic activity and not another gi¬ 
ant planet, which appears to be the most natural expla¬ 
nation. Hereafter we assume that the latter is true. From 
now on, we will also make the reasonable assumption that 
XO-2 was formed from a single, chemically homogeneous 
molecular cloud, which fragmented to form two proto- 
stellar nebulae, each with its own proto-planetary disk. 

The very strong correlation between abundance differ¬ 
ence and Tc shown in Figure complicates the simple 
explanation given above. Volatiles are just slightly defi¬ 
cient by about 0.015 dex while refractories are lower by 
about 0.090 dex in XO-2S compared to XO-2N. If metals 
were taken from XO-2S to make more planets, the effect 
is observed mainly in the refractories. Thus, based on 
our previous works, the planet population that explains 
this result cannot be the one detected by the transit and 
RV data. The planets detected by these traditional tech¬ 
niques have masses that are all greater than those ex¬ 
pected for rocky objects. In fact, they are all more mas¬ 
sive than Saturn. 

Even though the XO-2 planets likely have massive 
rocky, refractory-rich cores, they are also expected to 
be surrounded by large volatile rich envelopes. Based 
on the I6Cygni experiment, we expect the formation of 
these planets to produce mainly a constant offset in the 
abundance of ele ments with zero or just a very weak Tq 
correlation, as in Tucci Maia et al.| (2014). Thus, only 
the ~ 0.015 dex onset could be explained by the forma¬ 
tion of these planets, and not the Tq trend, or at least 
not most of it. 

Before discussing in more detail the effects described 
above, we must keep in mind that when interpreting 
chemical abundance anomalies we are using measure¬ 
ments of photospheric abundances, which are thought 
to be representative of the composition of the stars’ con¬ 
vective envelopes only. If the chemical signatures are 
imprinted when the stars’ radiative zones have fully de¬ 
veloped, i.e., after the first few million years since the 
stars’ birth, the amplitudes of these signatures are in¬ 
versely proportional to the mass of the stars’ convective 
envelopes. In other words, chemical anomalies due to 
accretion of dust-cleansed or metal-depleted gas will be 
diluted in stars with massive convective envelopes while 
they will be easily imprinted in stars with thin convec¬ 
tive envelopes. The XO-2 stars are somewhat cooler than 
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the Sun, with ~ 5500 K, but also significantly more 
metal-rich, at [Fe/H] ~ -|-0.4. For this combination of 
atmospheric parameters, which impl y masses very clos e 
to IMq, the evolutionary tracks by Siess et al. (2000), 
interpolated to the meta llicity of the system a rid usirig 
the solar abundances by Asplund et al. (2009) as refer¬ 


ence for the scaling, predict a convective envelope mass 
of ~ O.O3M0. 

5.2. Metallicity offset: gas giant planets 

Ignoring momentarily the additional depletion of re¬ 
fractories and concentrating only on the 0.015 dex over¬ 
all offset of metallicity between XO-2N and XO-2S, we 
can estimate the gas giant planet mass difference re¬ 
quired to produce it. Note that we refer to a mass differ¬ 
ence because both stars are known to host planets and 
thus metals have been taken away from each of them, 
but more so for the one that formed more giant plan¬ 
ets, i.e., XO-2S. Scaling to th e metallicity of the X O-2 
stars ([Fe/H] ~ 0.4), using the [Asplund et al.| (|2009D so¬ 
lar abundances as reference {(z/X)q = 0.0134, which 
leads to {Z/X)cz = 0.034 for XO-2), and adopting a 
convective envelope mass Mqz = 0.03 Mq for the XO-2 
stars, we can calculate the planet mass difference, pro¬ 
vided that we can estimate the me tallicities of the known 
planets, {Z/X)p. As discussed injRamfrez et al.| (|2011 
their Section 5.4), {Z/X)p is in the range trom 0.04 to 
0.12 for Jupiter. As a first order approximation, here we 
adopt [ZjXff = 0.1 for all XO-2 planets. Under these 
assumptions, we find that the mass difference that ex¬ 
plains a 0.015 dex offset in [X/H] is abou t 0.57 Mjup (we 
used t he correct form of Equation 1 in [Ramirez et al.| 
2011 ’ to figure out this number, Mp, which corresponds 
to the planet mass difference in this case). This value 
is well below the observed mass difference of the plan¬ 
ets detected by transit and RV data, which is at least 
1 Afjup (the latter value is the minimum difference be¬ 
tween the total mass of planets known around XO-2S, 
1.4Mjup -|-0.26Mjup, and XO-2N, O.OMjup). However, 
if we assume that (Z/X)p ~ 0.07, we can get the mass 
difference closer to the minimum expected value. 

The relative metal-depletion effect described above in¬ 
cludes the effect of dilution in the stars’ convective en¬ 
velopes. When we employed 0.03 Mq for the mass of 
these envelopes, we assumed that the planets were all 
formed after the stars had fully developed their radia¬ 
tive zones and had already thin convective envelopes. 
However, we know that planets, the gas giants in par¬ 
ticular, form within the first few million years from star 
birth, i.e., while the stars’ convective envelopes are more 
massive than seen during most of their main-sequence 
lifetimes. Another way of increasing the mass difference 
computed before to about 1 Mjup given the available data 
is by assuming that the planets formed when the stars’ 
convective envelopes were about 0.05 M©, roughly two- 
thirds of the way between the time the star was fully 
convective to the moment its convective envelope reached 
its final main-sequence mass of 0.03 M©. Given all these 
important considerations, the scenario described here to 
explain the 0.015 dex offset of volatiles and refractories 


for XO-2 is plausible 

5.3. Tc dependence: terrestrial planets 

The rocky bodies of the solar system (meteorites, as¬ 
teroids, terrestrial planets) are significantly enhanced 
in refractory elements (e.g., [McDonough fc Sun] 1995 


Alexander et al. [2001| [Zuber et al.j 2011)1 Interest- 
ingly, the anomalous chemical compositiori of the Sun 
relative to solar twins mirrors the abundance pattern of 
meteoritic-like and Earth-like composition rocks. In ad¬ 
dition, outside of the solar system, evidence of material 
with rocky composition is available from observations of 
white dwarfs po lluted by the re cent accretion of planetes- 
imals (see, e.g., [Xu et al.[[2014[ and references therein). 

In order to explain the upward Tc trend seen in Fig¬ 
ure [^ a depletion of or pollution by rocky, refractory- 
rich material is required. If indeed related to planetary 
material, the planets involved would have to be of the 
terrestrial class. To estimate the amount of rocky mate¬ 
rial needed to explain the abundance differences seen in 
Figure[^ we employed a dep letion model that follows the 
sche me of [Chambers [ ( 2010[ see also Melendez et al.|2012 
and [Ramirez et al.[[2U14a[) . We caicuiated the change 
in photospheric chemical composition that a star experi¬ 
ences due to the formation of rocks for a given mixture of 
Earth-like and meteoritic-like material, finding the best 
combination by cornparing the results of this calculation 
to the observationsWe find that about 20 M®, where 
is the mass of Earth, of an equal mixture of Earth¬ 
like material and meteoritic-like material is necessary to 
explain the data. This is a significantly larger amount of 
rocky material than the one required to explain the pecu¬ 
liar solar abundances (cs 4M m), for example (Melendez 
et al.[[2009 Chambers 2010). Nevertheless, this is natu¬ 
rally Scplained~5j7TIiem^er metallicity and the larger 
convective envelope mass of the XO-2 stars. 

Two different scenarios could explain the Tq correla¬ 
tion of the abundance differences seen between the XO-2 


® The correct version of that f ormula has (Mnv. -t- ly i-n) msrei 
of only Mqz in the denominator. |Rami'rez et al. pOlTI formula 
strictly valid only when Mp/Mci « 1. 


Mp) instead 
is 


stars. As in Melendez et al. (2009), we could interpret 
this abundance difference as a signature of the formation 
of 20 Mq of rocky material around XO-2S, the refractory- 
poor component. This interpretation does not necessar- 

t® The convective envelope masses that we are assuming are 
based on our estimates for the total masses of the XO-2 stars, 
which are both very close to 1 Mq . When estimating the masses 
of these stars, we employed our photometrically-derived Tgfr val¬ 
ues, which are about 140 K warmer than those measured by D15. 
If we use their atmospheric parameters instead, we would derive 
masses which are about 5 % lower, for which the convective en¬ 
velope masses are about 0.10 Mq. This would make it easier for 
us to explain the planet mass difference discussed here. However, 
we doubt that our stellar parameters are less accurate than those 
given by D15. 

Following|Chambers[ |[2010|l , for our calcul ations we adopted a 
mix o f the chemical composition of the Ea rth llWaenke fc Dreibus 
1988 1 and of CM chondritic meteorites flWasson &6 Kallemeyn 
19881. The rocky material was mixed into a convective enve¬ 
lope mass of 0.03 Mq, with photosp heric abundances obt ained 
by scaling the solar abundances from (Asplund et al.| (|2009|l . For 
each chemical element we computed the change m abundance as 
log[(Mcz + A4'rock,Barth+A4'rock,CM)/A^Cz], where Mcz is the mass 
of the convective zone, M^o^k Earth is the mass with Earth-like 
composition, and M^ock CM is the mass with CM chondritic com¬ 
position. We changed the Earth masses of both M^ock,Barth and 
.^rock CM until a satisfactory fit was found. Using only one type 
of material does not provide a good fit to both volatile and refrac¬ 
tory elements; both Earth-like and meteoritic-like composition are 
needed to match the observations. 
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ily imply that no rocky planets formed around XO-2N. 
However, if rocky planets did form around XO-2N, then 
the total mass of these planets would have to be lower 
by about 20 Mq than that of all rocky planets formed 
around XO-2S. In other words, more refractory-element 
mass needs to be taken away from XO-2S in order for 
this scenario to be qualitatively and quantitatively con¬ 
sistent with Figure Within this framework, we would 
need to explain why XO-2N formed no little planets, or 
fewer little planets than XO-2S. In any case, none of 
those planets have yet been detected. 

The second explanation that could be provided to ac¬ 
count for the Tc correlation is that previously-formed 
planetary material around XO-2N ended-up being ac¬ 
creted by the star at later times (not necessarily during 
the formation of star and planets). The amount of this 
material would have had to be also 20 , but if this also 

happened in XO-2S, then the amount of rocky material 
accreted by XO-2N would have had to be 20 M 0 larger 
than that for XO-2S. This idea was already proposed by 
D15 to explain the iron abundance difference that they 
found. In this case, one would need to explain why this 
late accretion of rocks happened only in one of the stars, 
or more so around one of the stars, namely XO-2N. Per¬ 
haps the fact that XO-2N has a gas giant planet in a very 
close orbit could be related to this phenomenon, since a 
migrating giant planet can help dragging super-Earths 
towards the host star. The transiting planet around XO- 
2N is indeed very close to its parent star, but note also 
that the two giant planets orbiting XO-2S are relatively 
close to their star, and have certainly migrated to their 
present-day orbits. 


5.4. Scenarios not involving planets 
A few studies have contradicted o ur previous obser¬ 


vational results (Schuler et al. 2011) or have provided 


we hf. 

Lve discovered Isee. e.g.. 

G 

onzalez Hernande 

et al. 

2010 

2013 

Onehag et al. 

201 

T 

Adibekyan et al. 

2014). 


Particularly relevant to this paper are two alternative in¬ 
terpretations, both unrelated to exoplanets. We address 

In Onehag et al. (2011), it is suggested that the pres¬ 
ence oFTiot''staJsin the immediate vicinity of the Sun 
and a solar twin in the M67 open cluster could have af¬ 
fected the final composition of these stars by radiative 
dust cleansing. While we cannot disprove this idea for 
the Sun and the solar twin in M67, it appears to be 
an unlikely scenario for XO-2. Even though the XO-2 
stars are widely separated (about 4 500AU), a massive 
hot star would have had to be much closer to only one 
of the two stars, XO-2S to be more specific, in order to 
radiatively eject dust from it while keeping the other star 
nearly untouched. Dust cleansing by a nearby hot star 
cannot easily explain refractory-element abundance dif¬ 
ferences for stars in binary systems and it probably does 
not work for the ca se of 16 Cygni eith er. 


In the work by Adibekyan et al. (2014) it is sug¬ 
gested that the birthplace of stars in the Galaxy may 
explain abundance peculiarities that correlate with the 
elements’ condensation temperature. They find that lo¬ 
cal old stars that have migrated from the inner disk have 
low refractory-to-volatile ratio. This is definitely not the 
reason behind our XO-2 result, where one of the compo¬ 


nents has a low refractory-to-volatile ratio relative to the 
other one. Under the very reasonable assumption that 
the XO-2 stars were born together and from the same 
gas cloud, they share a birthplace and are coeval, yet 
they exhibit an abundance anomaly of the same ampli¬ 
tude and of very similar nature as that seen in the Sun 
and a few other solar twins. This experiment as well as 
our previous work on 16 Cygni demonstrate that stellar 
age or birthplace cannot be the only parameters deter¬ 
mining the detailed chemical composition of stars and 
in particular the small anomalies that are detected only 
with ~ 0.01 dex precision abundances. 


5.5. Other binaries 

Binary systems like 16 Cygni and XO-2 are particularly 
interesting in the context of measuring highly-precise rel¬ 
ative abundances because each pair consists of two stars 
which are very similar to each other. Another pair dis¬ 
cussed in the literature in the context of p lanet signa¬ 
tures i s HD20782/81, which was analyzed by Mack et al. 
(2014). Planets have been detected around both compo- 


nents in th is binary. Co ntrary to our analyses of 16 Cygni 
and XO-2, |Mack et al.| found no difference in the chem¬ 
ical compositions of HD20782 and HD20781. We note, 
however, that the error bars of their relative abundances 
are relatively large (average A[X/H] = 0.04±0.07), most 
likely owing to the fact that the stars in this system are 
not twins of each other. Indeed, the Teff difference among 
them is greater than 400 K while in 16 Cygni and XO-2 
that difference is only about 60 K. 

In another recent study, our team has investig ated 
the chemic al composition of the HAT-P-1 system (Liu 


et al. 2014). There, the secondary star hosts a transit- 
ing planet of ~ 0.5Mjup, yet the two stars have iden¬ 
tical composition at the 0.01 dex level . One mig ht be 
tempted to use the HAT-P-1 result by |Liu et al.| as ev¬ 
idence against the hypothesis that gas giant planet for¬ 
mation imprints a signature in the chemical composition 
of the host star. However, we must note that as part 
of the hypothesis it is required that the stars’ convective 
envelopes are more massive than their final values during 
the formation of gas giant planets. Since gas giants form 
quickly, the dilution of their chemical signature by the 
more massive convective envelopes is facilitated. Thus, 
it could be that the planet host star in HAT-P-1 had a 
very massive convective envelope at the time its planet 
formed. This is not in contradiction with the 16 Cygni 
and XO-2 interpretations. In fact, it results naturally 
from the necessity to either have unusually long-lived 
disks or from the adoption of episodic accretion models 
of star formation, in which the specific episodic accretion 
history of each star determines the time evolution of the 
convective envelop e masses (see our ex tended discussion 
of these effects in Ramirez et al.|[20TT Section 5.4). 


5.6. Impact on Galactic Chemical Evolution 

Independently of their true nature, the elemental abun¬ 
dance differences detected between stars in binary sys¬ 
tems like 16 Cygni and XO-2 contradict one of the most 
important assumptions made in observational Galac¬ 
tic chemical evolution (GCE) studies, namely that the 
present-day composition of a star reflects that of their 
parent gas cloud. Gertainly, the magnitude of this effect 
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is not so large to introduce significant scatter into tra¬ 
ditional GCE-type plots of [X/Ee] versus [F e/H] such as 
those presen ted in the classic GCE paper by |Edvardsson 


et ah (1993). However, if 0.01 dex abundance analyses 


are done tor large samples of stellar twins, this process 
will become evident and must be taken into account. In 
particular, these effects may have a non-negligible impact 
on ongoing and future surveys designed to exploit “chem¬ 
ical tagging,” i.e., the association of groups of field stars 
according to their common compos ition, which would 


suggest a common site of formation (Freeman & Bland- 


ot our Galaxy. 


Hawthorn] 2002), as a tool for reconstructing the history 


6. FINAL REMARKS 

Stars in multiple systems are expected to have the same 
surface composition in the absence of accretion effects 
after the stars have developed thin convective envelopes. 
Thus, they remove the need to correct for Galactic chem¬ 
ical evolution effects which complicate the study of field 
stars. If, in addition, the two stars in a given binary are 
very similar to each other, the complications of standard 
chemical analysis become irrelevant if a careful strict dif¬ 
ferential approach is employed. Our careful spectroscopic 
analysis has revealed a clear difference in the chemical 
composition of the XO-2 binary system stars, for which 
a number of possible scenarios involving planets are pro¬ 
posed. 

The XO-2 experiment carried out in this paper pro¬ 
vides reliable observational information to further in¬ 
vestigate our hypotheses of chemical signatures due to 
planet formation as well as other alternative explana¬ 
tions. Naturally, high-quality observations of a larger 
number of twin-star binary systems, both for high- 
precision chemical analysis and exoplanet search, are 
highly encouraged to continue this research. At the same 
time, however, we should acknowledge the need to de¬ 
velop a population synthesis modeling approach to ex¬ 
plain all the available data in a more robust manner. It 
should be realized that observational results that char¬ 
acterize a single system cannot be used to fully sup¬ 
port or entirely reject our hypotheses. It is becoming 
increasingly clear that the interpretation of small chem¬ 
ical abundance anomalies in the context of exoplanets is 
much more complex than previously thought. Of course, 
this should come as no surprise. After all, we are using 
this information as indirect, left-over evidence of the very 
complex processes of star and planet formation. 
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